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The reflection and transmission amplitudes of waves in disordered multimode waveguides are stud- 
ied by means of numerical simulations based on the invariant embedding equations. In particular, 
we analyze the influence of surface-type disorder on the behavior of the ensemble average and fluc- 
tuations of the reflection and transmission coefficients, reflectance, transmittance, and conductance. 
Our results show anomalous effects stemming from the combination of mode dispersion and rough 
surface scattering: For a given waveguide length, the larger the mode transverse momentum is, 
the more strongly is the mode scattered. These effects manifest themselves in the mode selectivity 
of the transmission coefficients, anomalous backscattering enhancement, and speckle pattern both 
in reflection and transmission, reflectance and transmittance, and also in the conductance and its 
universal fluctuations. It is shown that, in contrast to volume impurities, surface scattering in quasi- 
one-dimensional structures (waveguides) gives rise to the coexistence of the ballistic, diffusive, and 
localized regimes within the same sample. 



I. INTRODUCTION 

The statistical properties of disordered systems is a 
rich and long-standing problem which attracts many ef- 
forts both theoretical and experimental. In quantum 
solid state physics much attention is paid to statistics of 
eigenfunctions and eigenvalues of closed disordered sys- 
tems. Disorder in solid state problems is usually rep- 
resented by impurities which are randomly distributed 
over the whole sample. For this sort of mesoscopic sam- 
ples with "bulk disorder" the number of well established 
statistical characteristics is enormous (see, for example, 
Rcf . [j] and references therein) . The success of mesoscop- 
ics may be ascribed to the existence of the non-linear o- 
model put forward by Efetova. Although it is a remark- 
able tool for studying mesoscopic effects, the u-model 
has, however, restricted validity. For example, the sys- 
tem size must be much greater than the mean free path. 
Generalization to chaotic ballistic systems (i.e. quantum 
billiards) has recently become a topic of a great interest. 
Progress in this direction has become possible due to tke 
recently proposed field theory for quantum ballisticsu. 
By exploiting the new tool, the authors of Refs. [| and 
^ managed to treat different correlators in a clean sys- 
tem within an extremely chaotic limit, when the typical 
relaxation time is of the order of the flight time (diffu- 



sive boundary scattering). They found that naive sub- 
stitution of the mean free path for the system size into 
correlators obtained for bulk disorder would give wrong 
results for the ballistic case, and that, in fact, systems 
with bulk and surface disorder are not equivalent. 

In this paper we address the problem of a ballistic sys- 
tem that is disordered in the sense that there are no bulk 
scatterers, and the only place where scattering occurs is 
at a rough boundary. We choose to explore not the statis- 
tics of the eigenvalues of completely closed systems (like 
resonators or quantum dots), but the statistical proper- 
ties of scattering and reflecting amplitudes in bounded, 
but open in one direction, systems (waveguiding struc- 
tures). The key entities for this problem are the trans- 
mission and reflection amplitudes. More specifically, we 
consider a TV-mode waveguide with the boundary corru- 
gated within a finite interval (L), and study the statis- 
tical properties of the transmission through, and reflec- 
tion from, the disordered segment of the waveguide. This 
problem arises naturally in the characterization of trans- 
port properties related to, for instance, optical waveg- 
uides and fibers, remote sensing, radio wave propagation, 
sonar, shallow water waves, and geophysical probingtirQ. 
On the other hand, it describes, .as well the electronic 
transport in mesoscopic systemaj'B, .Tpgdug especially rel- 
evant to conductance in nanowiresu'EJOtj. 
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The analogous problem with bulk scatterers was ad- 
dressed by many authors (see, for example, Refs. [| ||, 
|l4| and references therein). For a waveguide with bulk 
disorder all transmission coefficients T mn (sub-indexes n 
and m stand for number of the incident and transmit- 
ted mode, respectively) behave in a similar way. Due 
to the strong inter-mode mixing all information about 
the n-dependence is washed out after a few scattering 
events, which means that for L of the order of the scat- 
tering length i (and larger) all modes cross over from 
the ballistic regime to diffusion, and all channels become 
identical. As a result, there exists only one characteristic 
length scale for all transmittances, the so called local- 
ization length, that is believed to be equal to Ni. In 
the case L <C Ni, each T mn obeys Rayleigh statistics; 
as the length L increases, all channels (modes) undergo 
the same changes, and at L ^> Ni the crossover to the 
log-normal-di.st.ri b.n tinu (.typical for the localized regime) 
takes placey'^Olia-BeaH. It might seem that the only 
distinction of the problem with a rough surface from that 
with bulk disorder is that the scattering process takes 
place in a reduced effective volume which should lead just 
to a decrease of the mode mixing rate. Naive considera- 
tions would suggest that if we introduced a new localiza- 
tion length (which obviously must be much longer than 
that for the bulk scattering), all results well known for 
the "bulk" problem should be valid for the "surface" case 
after proper rescaling. However the situation is different 
and much more complicated. 

The goal of the present paper is to study the statisti- 
cal properties of waves transmitted through and reflected 
from a waveguide with rough boundaries. The length de- 
pendences of the reflection and transmission amplitudes 
for each realization of the surface profile are numeri- 
cally obtained by solving a system of linear differentis.1 
equations based on the invariant embedding equations^ 2 ]. 
Then we calculate the ensemble average and fluctuations 
of the reflection and transmission coefficients, reflectance, 
transmittance, and conductance. It is shown in this pa- 
per that the interplay between mode dispersion and sur- 
face scattering gives rise to many new and unusual (at 
least from the point of view of the intuition gained from 
studies of the volume scattering) effects; one of such ef- 
fects, the coexistence of different transport regimes at a 
certain length scale, has been previously reporteco. 

This paper is organized as follows. In Sec. || the the- 
oretical formulation leading to the invariant embedding 
equations for the matrices of the reflection and trans- 
mission amplitudes is developed. The description of the 
numerical implementation of those equations for the par- 
ticular waveguide geometry chosen here is detailed in 
Sec. [II. The results thus obtained for the average and 



fluctuations of the reflection and transmission coefficients 



are presented and discussed in Sec. IV, whereas those for 
the total reflection and transmission, and conductance, 
appear in Sec. The relevant conclusions derived from 



this work are summarized in Sec. VI 



II. THEORY 

We start from the wave equation 

(A + fc 2 ) W (R) = 0, (1) 

with the boundary conditions 

* (R = R s ) = 0, for x < and x > L, (2a) 
9* (R) 

= -^(R)-^^ i , forO<x<L, (2b) 

given on the unperturbed waveguide surface R — R s , 
which is translationally invariant along the x-axis (R = 
(x,r)). The boundary condition (||) corresponds either 
to a waveguide surface with a random admittance or 
to the Dirichlet boundary condition on a slightly per- 
turbed waveguide surface, £ denoting the random per- 
turbation. In the latter case the boundary condition (||) 
is an approximate one, containing the first two terms in 
the expansion of the exact (Dirichlet) boundary condi- 
tion about the unperturbed surface. 

Outside the region < x < L the solutions of the 
scattering problem under consideration have the form: 



(x, r) 
*™ {x,r) 




I Y (r)e~ ikmX t 

-; Am \ L J ° u r, 



x < 0, 



(3a) 



Xn (r) e 



-ik n x 



Xm (r) e 



x > L. (3b) 



The indexes "m, n" correspond to the outgoing and in- 
coming modes, respectively, and \n (r) are the eigenfunc- 
tions of the transverse wave equation 



(dr 2 



+ < Xn (r) = 0, 



(4) 



with k n = yjk 2 — k 2 j . 

By assuming that £n, the boundary condition 
can be also rewritten as 



*(R S ) + £(R S )$(R ; 
where the normal derivative 

"9* (R)" 



0, 



$ (Rs) = n (R^) 



<9R 



R,— R, s 



(5) 



9* 



is introduced, with n (R s ) being the normal to the un- 
perturbed surface R = R s . 

Let us employ Green's theorem in the form: 



*(R) 



3 D / 9 



d R 



<9R' 



* (R') dGo (R '' R) 



dR' 



9* (R') 

^^G (R',R) 



(7) 
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where the integral is taken over any region containing 
the point R and located inside the waveguide. The most 
convenient integration region is the unperturbed waveg- 
uide surface cut by two planes normal to the axis of the 
waveguide (they are included too). The left (right) cross- 
section must be placed to the left (right) of the point 
R = (x,r). The integral in Eq. (Q) can then be ex- 
pressed as the sum of the integrals over the cross-sections 
and the unperturbed waveguide surface with constraint 
< x' < L: 



* (R) = J dS' ■ W (R4) 



OGq (R',R) 



<9R 



(8) 



The integral I cr0 ss over the cross-sections may be easily 
calculated from Eqs. (||) and the Green's function 



Go (R, R') 



\ _ Xm (r) Xm (rQ c ik m \x-x'\ 



m—l 



2ik 



(9) 



<j> n (x, s) 



=n(r s ) 



dXn (r) 



dr 



The second result we derive from Eq. ( |TT| ) [by differ- 
entiating, setting R on the surface, and substituting the 
boundary condition (^)] is a closed equation for <& n (x, s): 

$„ (x, s) = cj> n (x,s) 



dx' I ds'i (x', s) $„ (x', s') G'q (x', x; s', s) . (14) 



Here G ' is the mixed normal derivative 

n " ( I I \ ® Go 1 V"> 1 dXndXn ik n \x-x'\ 

G (x ,x;s ,s = — — — = — > — — -e I ■ 

anon' 2i fc„ on an.' 

Differentiation of Eq. (Q) yields 

^ ^ 0m (x, s) CLran 



Then Eq. (§) takes the form: 



*„ (R) = -j=Xn (r) e 



-ik n x 



JdS'- *„ (H4) 



dG (R,R) 
9R 



(10) 



R'=R' 



The (oriented) surface element dS can be written as dS = 
ndS, dS = dx ds. The explicit form of the differential 
ds depends on the geometry under consideration (ds = 
rdcj) for circular cross-sections, dydz for rectangular cross- 
sections, etc.). Then we can rewrite Eq. (|To| ) as: 



* n (R) 



l 



L 



Xn (r) e 



-ik n x 



dx' <hds'* n {x ^ r , s) dGo(x>,x;r< s ,r) 



dn' 



(11) 



Two conclusions can be derived from Eq. (|TT|). First, the 
matrix of the reflection coefficients can be written [after 
substituting the explicit expression for Go from Eq. (|jj)] 
in the form: 

L 

Tmn = h J dx / ds ^ =( ^ m ^ e ^ ikmX ^ n rs ) ' ( 12 ) 



or [by the use of Eq. (ph] as: 



dx 1 & ds'i (x', s') d ^{x',s') G n ^ x _ ^ ^ 



2i 



dL 



ds'ct>* m (L,s)£(L,s)$ n (L,s). 



(16) 



By comparing Eq. (|l5|) with Eq. (|l4|) we obtain the rela- 
tion between the derivative <9<£> n (x, s)/SL and the func- 
tions {$ m (x, s)}: 



d$„ (x,s) 
dL 



$ m (x, s) a r , 



(17) 



With the aid of the latter equation, we now differentiate 
the matrix of reflection coefficients Eq. ( [l3| ) : 

dVrn.ri 



dL 



L< 

J dx <j> ds <j) m (x, s) £ (x, s) ^ $ m (x, s) (18) 



By substituting the explicit expressions for $ m (x, s) 
from Eqs. (|^b) and (^), and collecting all terms, we arrive 
at 



dr 



dL 2 
Here k = diag (k n ) and 



(19a) 



m = -i J dx I dscf> m (x, s) £ (x, s) (x, s) , (13) 
o 



where 



ds 4> m (s) £ (L, s) 0„ (s) 



0™ (s) 



n(r s ) 



dXn (r) 
9r 
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Analogous algebra leads to the equation for the matrix 
of transmission coefficients: 



dt i ^ , 
— = -te 
dL 2 



(19b) 



From the reflection and transmission amplitudes, we 
define the reflection and transmission intensities, respec- 
tively, 



(20) 



which yield the intensity coupled into the rath outgoing 
channel in reflection and transmission, respectively, for a 
given nth incoming channel. The reflectance and trans- 
mittance for the nth incident mode are 



(21) 



Finally, the total transmitted intensity in the case that 
all incoming channels arc incoherently populated, which 
is equivalent to the dimensionless conductance for elec- 
trons, is 



fJ 



(22) 



III. NUMERICAL CALCULATIONS 

For the numerical simulations we choose the simplest 
geometry (see Fig. |]): two parallel planes z = and z — d 
with ID deviations £ = £ (x) on one plane (z — 0) only, 
where £ is a ID stochastic process. Thus the transverse 
eigenfunctions acquire the form 



Xn (z) 



2 • ( \ - nn h 
— sin {K n z) , K n — " 



and the impurity matrix becomes 



,- 2 

71 • 



(23) 



(24) 



The 2N x 2N system of linear differential equations 
( |l9| ) is solved numerically by means of the 6th-order 
Runge-Kutta method. For each realization ^(x) (of 
length L max ) of an ensemble of randomly rough sur- 
face profiles, the matrices of reflection and transmission 
amplitudes are calculated as functions of the length L. 
These realizations obey Gaussian statistics (with 8 the 
rms height) with zero mean and a Gaussian correlation 
function 

W(\x - x'\) = 5~ 2 (£{x)Z(x')} = exp (~a- 2 (x - x') 2 ) , 

(25) 



where a is the transverse correlation length. The corre- 
sponding surface power spectrum is thus given by: 



g(Q)=^acxp(-(Qa) 2 /4) 



(26) 



The ensemble of surface realizations are numerically gen- 
erated as described in Ref. ^J. By averaging over N p 
such realizations, the mean values (A) and fluctuations 
5 A = ((A 2 ) — {A) 2 ) 1 / 2 of the relevant physical quantities 
are obtained. Hereafter we consider, unless otherwise 
stated, a waveguide of thickness d = 2.25A supporting 
N = 4 guided modes. 



IV. REFLECTION AND TRANSMISSION 
COEFFICIENTS 

In Fig. U the (T mn ) are shown (in a semi-log scale) for 
L < L max = 1500 A. Averaging was carried out over the 
results obtained for N p = 4000 realizations of the sur- 
face profile, whose roughness parameters are a = 0.2A 
and 5 = 0.03A. The asymmetry in the behavior of the 
different outgoing channels m is evident from this plot. 
The intensity of the incoming mode (T nn ) decreases with 
length for all n, this decrease being steeper the larger 
the transverse momentum k„, namely, the larger n. The 
transmission into other non-diagonal channels m =/= n 
also depends strongly on the mode m. In the beginning 
of the waveguide, this non-diagonal transmission slightly 
increases from zero, being stronger into higher modes m 
((T m+ i, n ) > (T mn ), with to+ 1, m ^ n). In this situation 
only single scattering is important, and we refer to this 
regime as quasi-ballistic (QB). In accordance with the re- 
sults of perturbation theory (PT) the intensity of mode m 
is proportional to the cross section for roughness- induced 
scattering from mode n into mode m, and to the length 
L\ as follows: 



(Trr 



2(5 f^n^m 



L 



d 2 (k nkm y/ 29ilkn - kml)L 

For the diagonal transmission, PT predicts 

(T nn ) = 1 - L^- ^Xn) 1/2 [9 ^ kn ~ k ^ 



(27) 



+g(\k n + k m \)} = 1 



L 

1QB- 



(28) 



In Fig. U the corresponding QB lengths l c ^ l from the 
preceding PT expressions have been plotted along with 
those obtained by fitting the numerical results shown in 
Fig. ^ to the expected linear functions, showing good 
agreement. 

The origin of such asymmetries lies in the surface-type 
disorder that randomizes the wave propagation through 
the waveguide. If we look at the impurity matri x (|24| ), 
which determines the scattering strength in Eqs. (|19[), it 
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is obvious that there are large quantitative differences in 
Vmn for distinct values of m and n. As a matter of fact, 
this matrix can be rewritten as: 



^^(fc„fc m ) 1/2 M m M„, 



where M n is 



(29) 



(30) 



Through a simple geometrical argument, as long as 5 <C 
d, M m can be considered the number of times that mode 
m hits (interacts with) the rough wall on its way along 
the waveguideEl In the case that d/X = 2.25, it turns 
out that, for instance, M4 w 8.5-Mi. This factor af- 
fects the impurity matrix not only for the outgoing mode 
through M m , but also for the incoming mode through 
M n . This gives a physically intuitive explanation of the 
results shown in Figs. |5[ and of all other processes that 
will be shown below. 

For larger L multiple scattering becomes relevant. This 
actually means that not only the scattering that brings 
energy to mode m from n should be taken into account, 
but also the leakage from m into other modes, as well 
as all interchanges between i and j for all As a 

result, the energy spreads over all modes: --diffusion (D) 
in the space of mode-numbers takes placec3. In fact, it 
is seen in Fig. || that all outgoing channels tend to yield 
comparable transmission intensities within the length of 
the plot, except for n = 1. 

Furthermore, after a long propagation distance 
through the waveguide so that mode conversion has suffi- 
ciently populated all outgoing channels, we observe that 
Tin > > T^n > T± n . This waveguide length is not 
reached within the length scale covered in Fig. g. Al- 
ternatively, we have increased the surface roughness to 
S = 0.1A, so that this result can be observed: the corre- 
sponding transmission coefficients, presented in Fig. ||, 
indeed confirm such behavior for a waveguide length 
L > 300A. Thus, the higher modes appear to be more 
strongly scattered. This is also manifested in the over- 
all behavior of the four outgoing channels depending on 
the incident channel. Figure || reveals that, beyond the 
waveguide length given above, the transmission curves 
appear to be qualitatively similar for all incoming modes 
n, but shifted downward as n is increased. The behav- 
ior of those transmission curves, following parallel expo- 
nential decays, is a signature of the onset of localization 
(L) due to the coherent interference of multiply scattered 
waves. (On the other hand, it is interesting to note that, 
if we zoom in Fig. |] for < L < 100A, the transmission 
curves are qualitatively similar to those of Fig. g.) 

Thus we have seen in Figs. B and || that the depen- 
dence of the impurity matrix ( p9|) on mode dispersion 
has significant quantitative consequences, and also strong 
qualitative consequences for the properties of wave prop- 
agation through surface-disordered waveguides. As has 



been demonstrated in Ref. £3|, it can give rise to an entan- 
gling of transport behaviors within the same waveguide 
length. In Fig. [|, the diagonal transmission coefficients 
T nn from Fig. ^ are shown in a log plot. The results have 
been fitted, where possible, to the well-known behaviors: 
QB as in Eq. (|28|); inverse power law expected for D 



l D 

frp \ ^ Tin . 



(31) 



and exponential decay associated with L. It is seen in 
Fig. I that, within the interval 10 < L/X < 70, QB 
transport of the (11) channel coexists with D for the 
(33) and (44) channels; also, D of the (11) mode coex- 
ists with L of the (44) channel for L/X f» 10 3 . This con- 
firms the coexistence of QB transport, D, and L predicted 
in Ref. ^ for 8-mode waveguides; nonetheless, in this 
4-mode waveguide the coexistence of all three regimes 
within the same length region is not observed due to the 
limitation in mode dispersion differences associated with 
the lower number of available modes. [On the other hand, 
it should be noted that our results, not shown here, re- 
veal such coexistence (QB-D or D-L) phenomena associ- 
ated with surface-type disorder in the case of narrower 
waveguides supporting only 3 or even 2 guided modes.] 
Interestingly, the impossibility of defining the D regime 
consistently for all outgoing modes at the same length 
scale makes irrelevant any comparison with theories such 
as the macroscopic, approach provided by random matrix 
theory (RMTjilH, which predicts (T mn ) = (g)/N 2 for 
all m, n. 

In addition, Figure [5] permits us to observe the 
crossover between different regimes for each mode sep- 
arately. For L < £4^, all 4 modes propagate almost 
ballistically. The transition from ballistic transport to 
D can be observed for all modes at the distinct waveg- 
uide lengths defined by the corresponding £® n (see Fig. Q 
below). Note that even though QB and D regimes ex- 
tend over different L regions, in both cases the regions 
are well defined by the magnitude of the transmission 
coefficient: (T nn ) w 1 for QB and (T nn ) m 10" 1 for D. 
This seems to indicate that, from the value of the av- 
erage transmission coefficient, the qualitative transport 
behavior can be roughly known, in agreement with Ref. 

although there exist remarkable differences concern- 
ing the length dependence and the entangling of regimes. 
Finally, coherent interference leads to L. In Figure [| all 
modes [mode (11) barely] reach the L regime within the 
maximum length of the waveguide L max . It should be re- 
marked that, whereas the exponential decay rate is sim- 
ilar for all n, the real onset of localization takes place at 
slightly different lengths: the lower n is, the longer the 
waveguide must be to observe L. 

The dependence of if^f , £® n , and £\\ n on surface rough- 
ness is shown in Fig. pi In this respect, with the aim of 
correctly defining l\\ n , the average of the logarithm of the 
transmission has been usee 
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(lnT nn )«-A (32) 

The predicted 8~ 2 behavior is seen in Fig. |(a) for the QB 
decay lengths, showing reasonable agreement with the 
PT results [cf. Eq. (p8|)]. Similar behavior is observed 
in Fig. ||(b) for £® n and £^ n . It is interesting to note 
that £® n is different for each n, whereas £% n coincides for 
all n. Thus the well known relationship £ L — N£ D is 
meaningless in this context. (Although £ L = N£y l seems 
to hold instead; in fact, it has been shown that if £R. is 
defined through the resistance, £ L = N£ D does holdEi) 

The normalized fluctuations ST rnn / (T mn ) are shown in 
Fig. (?]. It is evident that there are differences among 
the fluctuations for every channel, in agreement with the 
behavior of the mean values shown in Fig. ||; this corrob- 
orates the qualitative argument given above in connec- 
tion with the asymmetry in the mode scattering rates. 
Note that at the beginning of the waveguide, mode con- 
version into m ^ n leads to a variance of unity for the 
corresponding off-diagonal fluctuations, whereas the di- 
agonal m = n ballistic transport is revealed through the 
result that (6T nn ) 0. Furthermore, these diagonal fluc- 
tuations undergo the crossovers between QB, D, and L 
regimes as discussed above in light of the mean values 
(see Fig. ||). The (nn)-fluctuations exhibit an increase 
from towards 1 as the transport gradually changes from 
QB to D. The well known speckle pattern fluctuations 
5T nn / (T nn ) ss 1 for all (mn)-channels build up in the D 
regime, steadily increasing above 1 as the mode becomes 
localized. Therefore, the phenomenon of the QB-D and 
D-L coexistence can be recognized by comparing the di- 
agonal fluctuations with each other, corroborating the 
argument given above in light of the results for the mean 
values in Fig. For a sufficiently long waveguide, it can 
be seen in Fig. fj] that the normalized fluctuations tend 
to be larger the higher the outgoing-mode m is. A linear 
increase for all channels is observedEx For a given incom- 
ing mode n, the rate of increase is the same for all out- 
going channels m; nonetheless, the fluctuations appear 
to be larger the higher m is (within the noise accuracy). 
Analogously, the rate of increase is faster the higher n 
is. These considerations corroborate the arguments dis- 
cussed above on the mode selectivity of the scattering 
strength in connection with the transmission intensities 
in Fig. [|. 

Let us now turn to the study of the reflection coeffi- 
cients (R mn ). These are presented in Fig. |8| for the same 
waveguide considered in Fig. ||. The peculiar scattering 
properties induced by surface disorder manifest them- 
selves in an intricate manner in the reflection channels 
also. For sufficiently short waveguide lengths, we expect 
that the reflection coefficients should increase linearly as 
predicted by PT, through the expressions: 

{Rmn) - &lTS^ kn + kml)L 55 m (33) 

These PT QB decay lengths in reflection and those ob- 



tained from the numerical results are shown in Fig. |9|. 
The agreement is even better than in transmission, and 
the strong mode differences are indeed confirmed. Be- 
yond the QB regime for each incoming mode n, the di- 
agonal {R nn ) is enhanced as the waveguide length in- 
creases, whereas the remaining off-diagonal reflection 
coefficients exhibit differences with the following ten- 
dency: the higher is the mode m, the larger is the re- 
flection coefficient. Enhanced backscattering appears in 
the strong diffusive (or weak localization) regime as a 
result of the constructive interference of multiply scat- 
tered paths; nonetheless, the enhancement factor differs 
from 2 (predicted by simple arguments, provided that 
the single scattering contribution is absent or negligible). 
This comes as no surprise inasmuch as each channel may 
behave differently, as revealed in the transmission coeffi- 
cients (see Fig. through the entangling of transport 
regimes. Interestingly enough, this anomalous reflec- 
tion can result in an anomalous enhancement factored, 
larger than 2 [see Fig. ||(a): although the background 
cannot be unambiguously defined, (i?n) ~ 2.3(i? m i) for 
any m. I]. Therefore the reflection coefficients fail to 
satisfy^ (R mn ) = (I + S mn )N-^(l + N^^N - (g)). 
Figure || also seems to indicate that the onset of L does 
not introduce significant chan ges in the reflection coeffi- 
cients, in agreement with Ref. E9j. 

It is interesting to analyze the normalized fluctuations 
of the reflection coefficients (see Fig. |l^). Leaving aside 
the transient strong fluctuations for very short length 
scales (associated with the fact that the corresponding 
reflection coefficients are small), the diagonal m = n nor- 
malized fluctuations diminish with increasing length, this 
decrease being steeper the higher n is. Then they stabi- 
lize about the value 5R nn / {R nn ) = 0.5 as the D regime 
is reached, and remain constant when entering into the 
L regime. The off-diagonal normalized fluctuations [only 
the (42) channel is shown in Fig. |l(J), since all the rest 
are similar], on the other hand, remain about the vari- 
ance of unity linked to the speckle pattern fluctuations 
in reflection. Therefore weak localization halves speckle 
pattern fluctuations in backscattering. 

V. REFLECTANCE, TRANSMITTANCE, AND 
CONDUCTANCE 

We have thus seen that wave propagation along a 4- 
mode surface-disordered waveguide, due to the surface- 
type disorder, unlike for volume disorder, displays 
anomalous properties in the transmission and reflection 
coefficients as a consequence of the mixture of QB, D, 
and L regimes for different waveguide channels. Bearing 
in mind these properties, we now proceed to calculate 
the total transmission T„, reflection and dimension- 
less conductance g. 

In Fig. [n], we plot the mean total reflection {R n } 
and transmission (T„) coefficients, along with the trans- 
mission fluctuations 8T n , in our 4- mode waveguide with 
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8 = 0.03A. It is evident that these quantities differ sub- 
stantially from one incoming mode to another. The larger 
n is, the larger the mean reflectance and the smaller the 
mean transmittance (recall that energy conservation re- 
quires that R n + T n = 1). This could be qualitatively 
expected once again, at least in the limit of small waveg- 
uide len gth s, in light of the n-dependence of the impurity 
matrix (|29|), which is stronger for incoming modes with 
larger transverse momentum K n (higher n). 

The overall transport properties of the waveguide for 
a given incoming mode n can be understood through 
the behavior of the mean total transmission (T n ) [see 
Fig. |ll|(b)], as the summation of (T mn ) over all outgoing 
channels m. Figure [j^, which shows (T n ) in a log plot 
for a 4-mode waveguide analogous to that of Fig. |ll|(b) 
but with a rougher surface with 8 = O.IA, illustrates 
this discussion. To observe a definite transport regime 
in the total transmission, either the transport regimes 
of the different outgoing modes swap at certain length 
scales, or one of the (T mn )'s predominates over all oth- 
ers. Note that even though the most transparent mode 
gives the predominant contribution from a quantitative 
standpoint, it is not at all evident that the same is true 
for the qualitative behavior (for instance, a steeper, weak 
decay added to a larger, but smoother, background would 
yield as a result a quantity whose magnitude is of the or- 
der of the latter, but whose qualitative behavior is given 
by the former weak decay). In principle, it can be as- 
sessed that transport will obviously be QB for lengths 
shorter than ( name ly ; f or £, < £Q£, as Fig. 1^ re- 
veals through the QB linear decays (see also Fig. | ) ; this 
can also be verified in Fig. |ll|(b). Conversely, the expo- 
nential decay associated with L appears beyond waveg- 
uide lengths for which the lowest (In) mode is localized 
£f n , as seen in Fig. [TJ (see also Figs. |] and ||). Fi- 
nally, unlike the QB and L regimes, which must always 
be encountered for sufficiently small and large lengths, 
respectively, it is not obvious that the L~ l dependence 
(D regime) in the intermediate region between QB and 
L transport is observed. This effect is another manifes- 
tation of the entangling of different transport regimes of 
the (T mn ) due to the combination of surface-type disorder 
and large mode dispersion. In Fig. [ij, where the length 
dependence of (T„) is shown in a log scale, D should 
manifest itself through a linear decay [cf. Eq. (|3l])]. It 
is seen that this decay is practically absent for most in- 
coming modes. Only within a narrow waveguide length 
window for which the D length scales of the transmission 
coefficients swap, would the corresponding total trans- 
mission exhibit the expected L _1 behavior. In any case, 
it is obvious that the average reflection and transmis- 
sion coefficients fail to obey the predicted dependences 
(T„) = N-^g) and (-Rjirn N~ 1 (N - (g)) in the weak 
localization or D regimeoEj. 

All these transmission properties, stemming from the 
mixing of QB, D, and L transports produced by surface 
disorder, become even more pronounced in the dimen- 
sionless conductance g. Figure OI(a) shows a logarithmic 



plot of (g) for 4-mode waveguides with different surface 
roughness parameters S/X =0.02, 0.03, 0.04, 0.06, 0.08, 
and 0.1, whereas (lng) is plotted in Figure [l3](b) [and 
also ln(g}] for the two larger 8 values. The corresponding 
conductance fluctuations are given in Fig. [l4|. Follow- 
ing the argument mentioned above for the total trans- 
mission, now the QB regime is restricted to the shorter 
in this case £^ [see Fig. |l3](a)]. This is explic- 
itly shown in Fig. |^(a). Likewise, the true L behavior in 
the conductance is ensured for lengths beyond which the 
(11) channel appears localized; this is seen in Fig. |l3|(b) 
(in those cases for which £ L < L max ) through the lin- 
ear decay of (In g), and its departure from ln(g) (owing 
to the transition to log-normal statistics, for which the 
dominant contributions arise from the law-probability re- 
alizations that yield large conductanceso'EJ). Recall that, 
although the actual onset of L thus appears at slightly 
different lengths, the localization length for given rough- 
ness parameters as defined from Eq. ( |32"|) is the same for 
all modes (ran), and coincides with those for the trans- 
mittances and conductance [see Fig. |(b)] . On the other 
hand, the absence of the D regime in the conductance 
curves is explicitly displayed in Figure |l3|(a). Thus an 
anomalous conductance crossover from QB to L regimes 
is observed for the 4-mode, surface disordered waveguides 
with several roughness parameters used in obtaining the 
results of Fig. 13, giving additional confirmation of the 
predictions of Ref. |2^. Moreover, the conductance fluctu- 
ations within this anomalous transition regime stabilize 
in all cases shown in Fig. HI at a value (Sg « 0.29), 
which lies below the expected value of the quasi- ID uni- 
versal conductance fluctuations (UCF) for a well-defined 
D regime (Sg w 0.364, cf. Refs. fn| [TJ). This lower value 
of the UCF has been also numerically found in Ref. [ll], 
but no physical interpretation was given therein. When 
entering into the L regime, our results for the waveg- 
uides with rougher surfaces in Fig. [l4| reveal that the 
conductance-fluctuations decrease below the UCF region, 
as expectediJ'EJ. 

Two comments are in order concerning the anomalous 
QB-L crossover in the total transmission and conduc- 
tance mentioned above. First, it should be emphasized 
that surface disorder is not a sufficient condition. Ac- 
tually, in the case of surface-disordered waveguides with 
small mode dispersion and/or strong inter- mode mixing, 
so that the D-like regimes of different outgoing channels 
coexist, the L _1 diffusive dependence could also be ob- 
served. Nevertheless, even if such a D regime appears, 
our results still reveal an anomalous behavior, inasmuch 
as the mean total reflection and transmission fail to follow 
the predicted weak localization length dependences, as 
pointed out above. As a second remark, it is worth men- 
tioning that the D-like regime is enhanced in the average 
resistance (in which contributions from smaller transmis- 
sion coefficients predominate) , in contrast to the average 
conductance. The results presented in Ref. ^8] corrob- 
orate these comments, which thus show no discrepancy 
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with our results. 



VI. CONCLUSIONS 

The statistical transport properties of classical waves 
propagating along surface-disordered waveguides have 
been studied, with special emphasis on the distinctive 
imprint introduced by the surface-type disorder. For this 
purpose, the invariant embedding equations for the ma- 
trices of reflection and transmission amplitudes of the 
guided modes have been obtained. By means of Monte 
Carlo simulation calculations, in such a manner that for 
every surface realization the corresponding system of cou- 
pled differential equations is numerically solved, the sta- 
tistical quantities of interest are calculated. We have 
focused on the mean reflection and transmission coef- 
ficients, reflectances, transmittances, and conductance, 
along with their fluctuations. The interplay between 
mode conversion and the scattering processes produced 
by surface disorder results in new interesting effects in 
the physical quantities mentioned. 

For a sufficiently long waveguide, transmission into 
modes with lower transverse momentum (lower index m 
according to our notation) is favored (exhibiting smaller 
normalized fluctuations), no matter what the incoming 
mode n is. The influence of the incoming mode n is re- 
vealed in the decrease of all the transmission coefficients 
for higher n. For smaller waveguide lengths, we have 
confirmed, through the analysis of the length dependence 
of the mean transmission coefficients and fluctuations in 
the case of 4-mode waveguides, the entangling of ballis- 
tic, diffusive, and localized transmission of modes within 
the same waveguide region that was recently reported in 
Ref. ^3| in the case of 8-mode waveguides. 

With regard to the reflection coefficients, enhanced 
backscattering is observed when the rough waveguide is 
long enough, and the enhancement factor, as defined by 
the ratio mm[(R nn ) / (R mn )], can be larger than 2. In 
fact, the non-diagonal reflection coefficients tend to be 
smaller for the reflected modes with lower m for all in- 
coming modes n. The reflected speckle patterns exhibit 
reduced fluctuations in backscattering (SR nn )/(R nn ) < 
0.5, whereas the expected value of 1 is approximately 
obtained for other reflected channels (m 7^ n). Both av- 
erages and fluctuations behave similarly throughout the 
D and L regimes. 

The transmittance, namely, the normalized total en- 
ergy transmitted for a given incoming mode n, is larger 
for the lower modes n. It should be noted that, in spite 
of the small strength of the random component that is 
present on one of the waveguide planes, very strong re- 
flectances (of the order of or larger than 90%) can be ob- 
served for sufficiently long waveguides. This could be rel- 
evant in multimode, optical waveguides with spuriously 
rough boundaries over long propagation distances, where 
it constitutes an unwanted effect£9. 



We have also analyzed the effect of the entangling of 
QB, D, and L transport on the qualitative behavior of 
the mean transmittance and its fluctuations, showing an 
anomalous effective QB-L crossover. This has also been 
confirmed in the conductance calculations (average and 
UCF), for which the influence of the disorder strength 
has been shown. 

Finally, we would like to mention the very recent 
works by Garcfa-Martm et al on the diffusion-localization 
transitions and on the intensity distributionsEjtj, in 
nanowircs with surface-disordered hard-walls consisting 
of a number of slices with fixed length and random width, 
with similarities to the problem dealt with here. Their 
numerical results, based on a generalized scattering ma- 
trix formulation exploiting mode matching at each slice, 
exhibit also the non-isotropy of the scattering intensities, 
stressing however the agreement of the statistics at each 
transport regime with the RMT predictions. 

Experimentally, all these effects can be revealed in the 
transmission intensities through metal microwave guides, 
for which our theoretical boundary conditions apply very 
accurately. As pointed out in Ref. ^3], the appropriate ge- 
ometry would be a planar waveguide with two metallic 
plates, one of them at least randomly rough, with fea- 
sible dimensions and roughness parameters (as derived 
from the values used throughout this work upon scaling 
them by the wavelength in the centimeter range). Sim- 
ilar waveguides but with tube geometry have been suc- 
cessfully employed in connection with volume disordered. 
Also in the electromagnetic domain, optical waveguides 
or fibers (in the micron range) could be other experi- 
mental devicesEHI, accessible to such measurements, where 
the predictions of our calculations can manifest them- 
selves, although in order to make rigorous quantitative 
comparisons the boundary conditions might have to be 
revised. Furthermore, the propagation of acoustic waves 
or other classical waves through confined geometries with 
appropriate randomness can be adequately accounted for 
by our formulation, and thus similar phenomena might 
be expected therein. The conductance calculations can 
be also of interest in the electronic transport through 
nanowires. 
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FIG. 1. Illustration of the waveguide geometry. 
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FIG. 2. Mean transmission intensities (T mn ) as functions of 
length L in semi-log scale for a waveguide of width d/X — 2.25, 
supporting 4 modes, with disorder parameters a/A = 0.2 and 
S/X — 0.03: (a) incident mode n = 1; (b) n = 2; (c) n = 3; 
(d) n = 4. Averaged over N p = 4000 realizations. 
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FIG. 3. QB lengths l®^ (in wavelength units) in trans- 
mission versus outgoing channel m for the waveguide used in 
Fig. ^: Circles: n = 1; Squares: n = 2; Diamonds: n — 3; 
Triangles: n = 4. Filled (respectively, open) symbols denote 
the numerical simulation (respectively, perturbation theory) 
results. 
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FIG. 4. Same as in Fig. but for S/X = 0.1. 
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FIG. 5. Mean diagonal transmission intensities (T nn ) as in 
Fig. U in a log-log plot. Fits to linear (QB), L" 1 (D, with 
dots), and exponential (L) decays are shown. 
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FIG. 6. Typical decay lengths (in wavelength units) as 
functions of the height standard deviation <5 2 (in A 2 units), 
obtained from the mean transmission intensities (see text) 
for d/X — 2.25 and a/A = 0.2. Circles: n = 1; Squares: 
n = 2; Diamonds: n = 3; Triangles: n = 4. (a) i'JJf from 
numerical simulation data (symbols, asterisks denoting the 
conductance £® B ) and from perturbation theory (solid lines); 
(b) inn (filled symbols) and £^ n (open symbols), and £ L from 
the conductance (asterisks). 
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sities as functions of the length L, for the same parameters as 
in Fig. §. 
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FIG. 9. Same as Fig. ^ but in reflection. 
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FIG. 10. Same as Fig. |JJ but for the normalized fluctuations 
of the reflection intensities, including only the backscattered 
channels and the (42) off-diagonal channel. 
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FIG. 11. Mean total (a) reflection and (b) transmission 
intensities, and (c) normalized fluctuations of the total trans- 
mission intensities, as functions of the length L, for the same 
parameters as in Fig. B. 
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as in Fig. ]d, including fits to the QB and L regimes. 
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FIG. 13. (a) Mean dimensionless conductance as a func- 
tion of the length L in a log-log plot for a waveguide width 
d/X = 2.25, supporting 4 modes, with disorder parameters 
a/A = 0.2 and S/X = 0.02 (upper solid curve), 0.03 (dashed 
curve), 0.04 (long dashed curve), 0.06 (dot-dashed curve), 0.08 
(dotted curve), and 0.1 (lower solid curve). Averages over 
N p = 4000 realizations. Fits to the QB regimes are shown, 
(b) (In g) (without dots) and \n{g) (with dots) for S/X = 0.08 
(dotted curves) and 0.1 (solid curves), revealing the L regimes. 
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